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SUMMARY 


A  two-dimensional  finite-element  computational  procedure  is 
developed  for  the  accurate  analysis  of  the  strains  and  stresses  in 
adhesively  bonded  joints.  The  adhesive  layer  is  modeled  using 
Schapery's  nonlinear  single  integral  constitutive  law  for  uniaxial  and 
multiaxial  states  of  stress.  Effect  of  temperature  and  stress  level  on 
the  viscoelastic  response  is  taken  into  account  by  a  nonlinear  shift 
factor  definition.  Penetrant  sorption  is  accounted  for  by  a  nonlinear 
Fickean  diffusion  model  in  which  the  diffusion  coefficient  is  dependent 
on  the  penetrant  concentration  and  the  dilatational  strain.  A  delayed 
failure  criterion  based  on  the  Reiner-Wei senberg  failure  theory  has  also 
been  implemented  in  the  finite  element  code.  The  applicability  of  the 
proposed  models  is  demonstrated  by  several  numerical  examples. 
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1.  INTRODUCTION 


1.1  Preliminary  Comments 

During  the  last  two  years  of  research  for  the  office  of  Naval 
Research  a  finite-element  computer  program,  NOVA,  has  been  under 
continuous  development  for  an  accurate  analysis  of  adhesively  bonded 
joints  (see  [1,2]).  The  finite  element  program  NOVA  uses  linearly 
elastic  elements  to  model  the  adherends.  The  adherends  can  be 
represented  as  isotropic,  orthotropic  or  laminated  composite 
materials.  The  large  displacements  and  rotations  experienced  by  the 
adherends  in  many  types  of  loading  are  accounted  for  by  the  updated 
Lagrangian  description  of  motion.  The  adhesive  layer  is  modeled  using  a 
special  element  that  employs  a  multi-axial  extension  of  Schapery's 
nonlinear  single  integral  stress-strain  law  as  the  constitutive 
equation.  The  finite  element  formulation  based  on  the  updated 
Lagrangian  incremental  equations  is  presented  in  detail  in  (1,21.  The 
element  library  contains  an  eight  noded  isoparametric  element  which 
employs  quadratic  interpolation  functions  to  represent  the  displacement 
field  as  well  as  element  geometry.  The  program  can  be  used  to  conduct 
plane  stress,  plane  strain,  or  axisymmetric  analysis  of  an  adhesively 
bonded  structure  subject  to  a  time  varying  thermal  and  mechanical 
loading.  The  program  is  modified  recently  to  include  a  nonlinear 
Fickean  moisture  diffusion  model  and  a  energy-based  delayed  failure 
criterion.  This  report  describes  the  theory,  finite  element  formulation 
and  applications  of  NOVA  to  moisture  diffusion  and  delayed  failure 
problems. 
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1.2.  A  Review  of  the  Literature 

Over  the  years  several  time-dependent  failure  criteria  have  been 
proposed  for  predicting  yield  and  fracture  of  polymeric  materials. 

Nagdhi  and  Murch  [3]  and  Crochet  [4]  have  used  a  modified  von  Mises 
criteria  for  viscoelastic  materials  by  assuming  that  the  radius  of  the 
yield  surface  depends  upon  the  strain  history.  An  energy  based  delayed 
failure  criterion  for  polymeric  materials  subjected  to  isothermal  creep 
was  developed  by  Reiner  and  Weissenberg  [51.  According  to  this  theory, 
failure  occurs  when  the  stored  deviatoric  strain  energy  in  a 
viscoelastic  material  reaches  a  certain  maximum  value  called  the 
resilience,  which  is  a  material  constant.  Bruller  [6]  and  Hiel  et  al. 
[7|  applied  the  Reiner-Weisenberg  failure  criterion  to  various 
viscoelastic  materials,  including  composite  laminates,  and  obtained  good 
agreement  with  experimental  observations. 

It  is  now  well  known  that  moisture  diffusion  can  have  a  significant 
effect  on  the  stress  field  within  an  adhesive  layer  in  a  bond.  Weitsman 
(8)  used  a  variational  method  coupled  with  Fickean  diffusion  law  to 
study  the  interfacial  stresses  in  viscoelastic  adhesive-layers  due  to 
moisture  sorption.  From  the  results  of  this  uncoupled  linear 
thermoelastic  analysis,  he  concluded  that  the  location  of  the  maximum 
interfacial  tensile  stress  depends  on  the  geometry  of  the  joint  as  well 
as  the  progress  of  the  diffusion  process  within  the  joint.  Weitsman  [9] 
used  the  correspondence  principle  to  generate  a  linear  viscoelastic 
solution  from  the  linear  elastic  analysis  of  moisture  sorption  within  an 
adhesive  layer.  He  observed  that  the  viscoelastic  analysis  predicts 
detrimental  effects  due  to  stress  reversals  caused  by  fluctuations  in 
relative  humidity,  that  are  not  predicted  in  an  elastic  analysis. 
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However,  he  acknowledged  the  omission  of  the  effect  of  moisture  content 
on  the  viscoelastic  response  of  the  resins  in  his  analysis. 

Tobing,  et  al.  (10]  used  the  finite  element  method  to  study  the 
micro-mechanical  effect  of  moisture  sorption  in  graphite-epoxy 
composites.  Using  a  constitutive  equation  based  on  the  Flory-Huggins 
lattice  model  for  polymer  solvent  interactions,  they  concluded  that  the 
stresses  at  the  graphite-epoxy  interface  have  a  strong  dependence  on 
moisture  content,  fiber  spacing,  and  applied  load. 

Yaniv  and  Ishai  [111  developed  a  linear  viscoelastic  closed  form 
solution  as  well  as  a  nonlinear  finite  element  solution  algorithm  to 
study  the  hygrothermal  effects  in  a  bonded  fiber-reinforced 
plastic/aluminum  system.  The  linear  solution  was  used  for  short-term 
predictions  at  low  strain  levels,  whereas  the  finite  element  solution 
was  used  for  long  term  predictions  in  which  geometric  and  material 
nonlinearities  were  taken  into  account.  The  authors  observed  that  the 
presence  of  moisture  tends  to  considerably  reduce  the  stress  level  in 
the  adhesive  layer  and  may  lead  to  significant  variation  in  the  time- 
dependent  deformation  of  the  test  specimen  as  compared  to  the  dry  state. 

In  the  references  cited  above,  various  authors  have  underscored  the 
effect  of  moisture  content  on  the  viscoelastic  response  of  the  test 
specimen.  However,  the  effect  of  the  viscoelastic  stress  field  on  the 
diffusion  coefficient  was  not  considered.  Lefebvre  et  al.  [12]  extended 
the  free  volume  concept  to  define  a  diffusion  coefficient  that  is  a 
function  of  temperature,  dilatational  strain  and  solvent 
concentration.  The  proposed  nonlinear  diffusion  model  showed  good 
predictive  capability  for  different  values  of  temperature  and  moisture 
conentrations.  They  concluded  that  in  order  to  obtain  an  accurate 
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solution  for  the  hygrothermal  effects  within  an  adhesive  bond,  the 
nonlinear  diffusion  problem  needs  to  be  solved  in  conjunction  with  the 
nonlinear  viscoelastic  boundary-value  problem  until  convergence  is 
achieved. 


2.  NONLINEAR  VISCOELASTIC  FORMULATION 

2.1  Introduction 

A  thermodynamically  consistent  theory  for  a  single  integral 
representation  of  nonlinear  viscoelasticity  was  first  proposed  by 
Schapery  [131.  The  law  can  be  derived  from  fundamental  principles  using 
the  concepts  of  irreversible  thermodynamics.  A  comprehensive  review  of 
the  thermodynamics  basis  of  Schapery 's  theory  has  been  presented  by  Hiel 
et  al .  1 14] . 

The  following  two  sections  deal  with  the  review  and  application  of 
Schapery 's  single  integral  constitutive  law  to  problems  with  uniaxial 
and  multiaxial  states  of  stress,  respectively.  The  constitutive 
equations  thus  obtained  are  suitable  for  non-linear  viscoelastic  finite 
element  analysis. 

2.2  Uniaxial  Stress  State 

The  uniaxial  nonlinear  viscoelastic  constitutive  equation  of 
Schapery  (13|  can  be  written  for  an  isotropic  material  as, 

^  gj  /  a0(/  .  ,5^  ^  jg^.^lds.  (2.1) 

In  Eq.  (2.1),  represents  uniaxial  kinematic  strain  at  current  time 

t,  is  the  Cauchy  stress  at  time  t,  0^  is  the  instantaneous  elasti' 

compliance  and  aO(iii)  is  a  transient  creep  compliance  function. 

Superscript,  t,  denotes  current  time.  The  factor  q^  defines  stress  and 

^0 
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temperature  effects  on  the  instantaneous  elastic  compliance  and  is  a 
measure  of  state  dependent  reduction  (or  increase)  in  stiffness, 

^0  ^  T'^^f’sient  (or  creep)  compliance  factor  has  similar 

meaning,  operating  on  the  creep  compliance  component.  The  factor 
g^  accounts  for  the  influence  of  load  rate  on  creep,  and  depends  on 
stress  and  temperature.  The  function  ii>^  represents  a  reduced  time  scale 
parameter  defined  by. 


=  J  (<t) 


'ds. 


(2.2) 


where  a^-j.  is  a  time  scale  'shift  factor'.  For  thermorheologically 
simple  materials,  a  =  a(T)  is  a  function  of  temperature  T  only.  This 
function  modifies,  in  general,  viscoelastic  response  as  a  function  of 
temperature  and  stress.  Mathematically,  a^^.  shifts  the  creep  data 
parallel  to  the  time  axis  relative  to  a  master  curve  for  creep  strain 
versus  time.  In  this  model,  four  material  parameters  (g^,  g^,  g^  and  a^) 
are  available  to  characterize  the  nonlinear  behavior  instead  of  only  one 
with  the  time  scale  shifting  procedure  of  Knauss  and  Emri  (15). 

The  transient  creep  compliance,  AD(i')»  can  be  expressed  in 
exponential  form  as, 

-X 

aD(u.)  =  I  D^(l  -  e  I,  (2.3) 

r 

where  0^  and  x^  are  constants.  Substitution  of  Eq.  (2.3)  in  Eq.  (2.1) 
gives. 


t  s 

=  gX"  ^  ^5  J  I  0,11  -  e  ^  Ig^'lds. 

0  r 


(2.4) 


Letting  the  product  g^o^  be  expressed  as  6^  and  simplifying  the 
integrand  on  the  right  hand  side  of  Eq.  (2.4)  yields. 
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g^D  +  gf  y  D  J  ^  G^ds  -  gf  T  D  /  e  ^  ds. 

^0  0  ^1  ^  r  •'q  ds  1  ^  r  •'q  ds 


(2.5) 

The  third  integration  term  on  the  right  hand  side  of  Eg.  (2.5)  is 
now  separated  into  two  parts,  the  first  part  having  limits  from  zero 
to  (t  -  At)  and  the  second  integral  spanning  only  the  current  load  step, 


i.e.,  from  (t  -  At)  to  t.  Hence, 


t  -  n,^)  ^gS 

J  ®  h7- 


t-At  -X^(il.^-4,^)  ^gS 


t  -X  pS 

+  J  e  ^  ^  ds.  (2.6) 

t-At 


The  first  term  on  the  right  hand  side  of  Eq.  (2.6)  can  be  rewritten  as. 


e  35-  ds  =  e  , 


(2.7) 


where 


.  t  t  t-At 
Au)  =  -  ii) 


(2.8) 


(2.9) 


The  second  integral  on  the  right  hand  side  of  Eq.  (2.6)  is  now 
integrated  by  parts.  To  carry  out  the  integration,  it  is  assumed 
that  G^  varies  linearly  over  the  current  time  step  At.  Hence, 


dG^  .. 

J  6  —  ds 

t-At 


-X  (ili''-il>^)  ^ 

dG^  e 

ds  x^  * 


t  ,2„s 

/  -  ds 

t-At  ds'^  V 
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1. 


(2.10) 


dG'^  ,1 


In  arriving  at  the  second  step,  it  is  assumed  that  is  linear  and 
hence  its  second  derivative  is  zero.  Since  G^  has  been  assumed  to  be  a 


linear  function  of  time  over  the  current  load  step,  we  can  write 

dG^  _  G^  -  G^‘^^ 


(2.11) 


Substitution  of  Eq.  (2. 11)  into  Eq.  (2.10)  gives 


-K  )  (JqS  .  t 

e  ^  ds  =  [G  -  G  |B^, 


(2.12) 


where. 


„t  I  -  e‘V^'^ 
^  =  - 1 - 


(2.13) 


Substituting  Eq.  (2.9)  and  Eq.  (2.12)  back  into  Eq.  (2.5),  and 
writing  obtains 


■  I5X  •  5^)  i  “r  - 


^t  V  n  In  t-At^t 

q ,  ;  D  q  T  3  j 
^  1  /  r '  ^  2  r 


t-At^t  t-f,t  ' 

3^j  -  e  q^  If. 


(2.14) 


Defining  instantaneous  compliance  Dj  as  the  compliance  term 
multiplying  the  instantaneous  stress  o'",  and  the  remaining  terms  in  Eq. 
(2.14)  as  hereditary  strains  E*",  one  can  write 


=  q\<3^  +  E^  :  F(o), 


(2.15) 


where 


"l  '  9^  *  9I9)  1  “r  -  9(9)  I  V) 
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e’'  =  9\{l 


(2.17) 


Equation  (2.15)  expresses  Schapery's  single  integral  constitutive 
law  in  terms  of  a  stress  operator  that  includes  instantaneous  compliance 
and  hereditary  strains. 

It  is  to  be  noted  that  the  term  in  Eq.  (2.17)  is  the  r^^ 

component  of  the  hereditary  integral  series  at  the  end  of  the  previous 
load  step  (i.e.  at  time  equals  t  -  At).  The  expression  for  the 
hereditary  integral  at  the  end  of  the  current  load  step  (i.e.  at  time  t) 
can  be  derived  in  the  form  of  a  recurrence  formula.  From  fiq.  (2.9),  one 
can  write 

'  sr 


t-At  .pS  t  .pS 

=  J  e  ^  i-ds  -  J  e  ^ds.  (2.18) 

0  t-At 


Using  the  results  from  Eqs.  (2.9)  and  (2.12),  the  above  equation  can  be 
reduced  to 


t  t'At  ^  i„t  t  „  t-At  t-At,„t 

^  -  92 


o  l6^. 


(2.19) 


where  6^  is  defined  by  Eq.  (2,13). 


2.3  Multj-dxial  Stress  State 

For  a  thermorheological ly  simple  anisotropic  viscoelastic  material 
under  a  multiaxial  state  of  stress,  the  constitutive  law  proposed  by 
Schapery  1131  is. 


3G, 


^ij  ■  ~  3a 


IJ 


3o 

-  Ae 

3o..  mn 


(2.20) 
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Ae  =  f  -  idM  —  (a.  ./a.)d 

mn  •'  mn'  '  3t  '  ij 


(2.21) 


where  e^j  and  a^j  are  the  components  of  the  strain  and  stress  tensors, 

respectively;  Gn  is  the  Gibbs  free  energy,  o.  .  and  are  the 

1  j  mn 

components  of  the  second  and  fourth  order  material  property  tensors, 
respectively;  and  ag  is  a  material  kernel  function  defined  in  (13).  The 
quantities  Gp,  ag  and  are,  in  general,  functions  of  ten 
variables,  a^.^.  and  temperature  f.  Note  that  all  repeated  indices  in 
Eqs.  (2.20)  and  (2.21)  are  to  be  summed  over  their  range  (1,2,3). 

Due  to  the  complex  nature  of  Eq.  (2.20)  it  is  not  possible  to 
determine  the  material  properties  in  this  equation  from  the  uniaxial 
tests  outlined  in  |131.  However,  it  is  possible  to  construct  a  set  of 
small  strain,  three-dimensional  constitutive  equations  from  Eq.  (2.20), 
which  is  consistent  with  the  thermodynamic  theory  in  [131  and  yet 
enables  all  properties  to  be  evaluated  from  uniaxial  tests.  The 
assumptions  which  need  to  be  made  for  this  purpose  are  as  follows: 

(a)  The  Gibb's  free  energy  Gp  is  a  quadratic  function  of  stress, 

(b)  a.j  =  ...  (2.22) 

When  the  free  energy  Gp  is  a  quadratic  function  of  stress  one  has 


— ^  =  S"""(0)o  , 

Ij'  '  mn’ 


(2.23) 


where  3.j(0)  are  the  instantaneous  components  of  the  linear  viscoelastic 
creep  compliance  tensor.  Equation  (2.23)  imp-lies  that  the  initial 
response  of  the  material  is  linearly  elastic  under  suddenly  applied 
stresses,  which  is  often  the  case  for  metals  and  plastics. 

The  second  assumption,  on  the  other  hand,  leads  to  the 
linearization  of  the  coefficient  of  the  transient  term  in  Eq.  (2.20). 
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Mathematically,  this  is  given  by, 


3o  1,  if  i  =  m  and  j  =  n 

_ mn  ^  j 

^°ij  0,  if  i  *  m  and  j  *  n 


(2.24) 


Equation  (2.24)  implies  that  the  jump  in  strain  due  to  load  application 
equals  the  jump  when  the  load  is  removed.  This  behavior  is  exhibited  by 
some  types  of  plastics  (16). 


Substituting  Eqs.  (2.23)  and  (2.24)  in  Eqs.  (2.20)  and  (2.21),  one 


obtains 


“(j  -  /  -  *’)  It 

—GO 

where  g^  =  Equation  (2.25)  is  a  set  of  three-dimensional 

constitutive  equations  for  anisostropic  viscoelastic  materials  which 
includes  the  nonlinear  functions  g^  and  a^y  appearing  in  the  uniaxial 
relations  (2.1)  and  (2.2).  Note  that  the  functions  g^  and  a|^^  are 
expressed  as  a  function  of  the  octrahedral  shear  stress. 

For  a  homogeneous  isotropic  nonlinear  viscoelastic  material,  Eq. 


(2.24)  reduces  to  the  form  presented  by  Schapery  ([13]  and  1161), 


Aj  =  ^  'D  -  JWg 


2"mm  ' i j’ 


(2.27) 


where 


iJl.g^oJj)  =  J(0)o.j  +  -  /)  (gJc;.)dT 

ID  -  -  (0(0)  -  J(0)laj^^  ^  [aD(i^  -  ;^) 


(2.28) 


-  T7 


(2.29) 
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(2.30) 


^  t  , 

4,  -  4i'  =  f  (a^-p)  ds  . 

Expanding  Eq.  (2.27)  term  by  term  for  the  strains,  one  has 

efl  =  {D}{g2'’5l^  ■  ‘^■^5242^  •  (2.31a) 

Simi larly , 


622  =  (D  -  ^  ':o}tg^22^ 

y\^  =  2iJHg^c;2) 

e^3  =  tO  -  JHg^^i'i  ^  (D)  (9^0^31.  (2.31b) 

The  transient  components  of  the  creep  and  shear  compliances  can  be 
written  in  the  form  of  the  Prony  series  as, 


:0(.)  =  V  D^ll  -  e  '^*1 

r 

aJ(.)  =  ;  J^ll  -  e  '■'l,  (2.32) 

r 


where  and  are  the  reciprocal  of  the  retardation  times  in  creep  and 
shear  respectively.  Also,  let. 


0(0)  = 

J(0)  =  (2.33) 

Considering  a  term  of  the  form  D-g^^^jl  in  Eq.  (2.31)  and  substituting 
Eqs.  (2.32)  into  Eq.  (2.33)  gives. 


(Dfig^^jj 


=  Do4j  *  4  I  - 


192'"] 


(2.34) 


Recognizing  that  Eq.  (2.34)  is  similar  in  appearance  to  Eq.  (2.6)  and 
making  use  of  the  results  derived  in  Section  2.2,  one  obtains 
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(2.35) 


where  Dj  is  the  instantaneous  creep  compliance  function  at  time  t, 

-  ej). 

Q^j  are  the  hereditary  strain  components  due  to  tensile  creep  at  time 


AiJ) 


ot  V  I„t-At„t  t-At  „  r  „t-At, 

«lj  •  I  192  ViJ  -  ' 


-X  Ai/ 

8^  -  --  g  - 

X  A.' 

r 

and,  .  .  are  components  of  the  hereditary  integral  given  by  the 
^ '  J 

recurrence  formula, 

,  ,  t 

^r,ij  ^r,ij  ‘^2  1 j  ^2  ij  '  r  • 

Similarly,  a  term  of  the  form  'J;{g^c^j;  in  Eq.  (2.31a)  can  be 
expressed  as. 


(2.36) 


where,  is  the  instantaneous  shea’'  compliance  function  at  time  t. 


j)  •  Jo  ♦  9)  I  -V"  -  V>’ 


are  the  hereditary  strain  components  due  to  shear  at  time  t, 

,  t 

pt  V  .  t-At,t  t-At  nt-fit, 

r..=)q„  o,.  -e  p  .. 

U  ^  ^2  r  ij  ^r,ij' 


-t  1  -  e 
■  r  .  t 
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and,  .  .  are  components  of  the  hereditary  integral  given  by  the 
^  ^  J 

recurrence  formula 


-"r**  t-At  .  I  t  t 
'  Pr.lJ  "  'Vlj 


-  g 


t-At  t-Ati ^t 
2  °ij  '^r 


Substituting  Eqs.  (2.35)  and  (2.36)  in  Eqs.  (2.31)  and  dropping 
superscripts,  one  obtains 


^  "  '^I^“22  ■  '^I^°33  ^  ^11 

®22  ^  ^1°22  ■  '^I^°33  ^22 

M2  =  ^'^I'’12  ^  ^12 

®33  ~  (^j  ~  ~  '^I^°22  ^  ^l'^33  ^  ^33*  (2.37) 


where 


^1 

=  ^11  ^  ^22  ^  *^33  ■ 

p 

22 

■  '’33 

^22 

=  ^11  ^22  ■"  *^33  ■ 

'’ll 

'’33 

^12  "  ^''12 

^33 

=  +  Q22  *  Q33  - 

'’ll 

'’22* 

(2.38) 

Equation  {2.37)  can  be  e^Pi^^ssed  in  matrix  form  as 

{e^  =  1N1{t  +  ;H;.  (2.39) 


Note  that  the  left  hand  side  of  Eq.  (2.39)  is  a  vector  containing 
the  algebraic  difference  of  kinematic  strains  (e}  and  dilatational 
strains  fs^j9). 
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{(^11  ■  9  )  »  Y  ^  5  I  (  t -j-}  ~  9)/, 


'22 


12’^''33 


(2.40) 


while  (o)  contains  four  components  of  Cauchy  stress, 


(o}  =  1 1  ’ ^22 ’ "^12 * *^33 ^  ’ 


and  (H}  is  a  vector  of  hereditary  strains  given  by, 

(H}  = 


The  matrix  '^'l  is  a  4  x  4  coefficient  matrix  given  by. 


(N|  = 


Pre-multiplying  Eq. 


“I 

(Dj-Jj) 

0 

(Dj-Jf) 


(Dj-Jj)  0  (Dj-Jj) 

Dj  0  (Dj-Jj) 

0  2J^  0 

(Dj-Jj)  0  Dj 


(2.39)  by  an  explic't  expression  for  stresses 


in  terms  of  strains  is  obtained: 

=  lMl(  e-  -  H;)  ,  [M|  =  !Nl'\ 


(^.-ll) 


Equation  (2.41)  provides  a  general  viscoelastic  constitutive 
relation  that  can  be  applied  to  either  plane  stress,  plane  strain  or 
axisymmetric  problems.  For  plane  strain,  the  out-of-plane  strain 
component  e^^  is  identically  zero.  The  corresponding  stress 
component,  may  be  obtained  from  Eq.  (2.41)  by  setting  633  =  0. 

Since  for  the  plane  stress  case  is  identically  zero,  the 

corresponding  strain  component  633  can  be  evaluated  from  Eq.  (2.39)  as 

'33  ■  <”1  -  *  ’??>  •  "33- 

Note  that  the  use  of  creep  and  shear  compliances  as  material  property 
input  allows  the  Poisson's  ratio  to  be  time-dependent.  Hence,  the 
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prsent  formulation  is  applicable  to  any  thermorheologically  simple 
isotropic  viscoelastic  material  over  any  length  of  time. 

For  the  special  case  where  the  Poisson's  ratio  is  a  constant  with 
time,  then 

=  (1  +  v)D(4^).  (2.42) 

The  matrix  [N]  takes  the  form. 


and,  the  corresponding  hereditary  strains  are, 

Hii  =  ^  *^33^ 

^22  "  ^22  ■  -^^11  ^ 

^2  ^  ^  ■■)'^12 

H33  =  Q33  -  -(Qii  ^  Q22)  (^•'’3) 

If  the  viscoelastic  properties  of  a  material  are  defined  by  its 
bulk  and  shear  compliances  instead  of  the  creep  and  shear  compliances, 
then  the  creep  compliance  0(.)  in  Eg.  {2.27)  is  replaced  by  the  bulk  and 
shear  properties.  Using  the  viscoelastic  relationship  betwen  creep, 
bulk  and  shear  compliances  given  by, 

DCv)  =  5  M(0  +  I  J(«.),  (2.44) 

and  substituting  in  Eg.  (2.27),  the  relation  in  Eg.  (2.39)  is 
obtained.  However,  for  this  case  the  matrix  |N|  has  the  form. 


15 


[N] 


<5"l  *§Jl> 

(5"l  - 


(Am 


0 

I  ■  3 


T  Jl) 


(i«,  -ij,)  0 

0 

0  2J^ 

(^Mj-iji)  0 


(g  Mj  -  3  Jj) 

^9  "  i 

0 

(g  Mj  +  f  Ji)j 


(2.45) 


and  the  corresponding  hereditary  strains  are, 

^11  "  i  ■"  ^22  ^33^  I  ^11  ■  I  ^'^22  ^33^ 

^22  "  9  ^22  ^33^  3  ^^22  '  I  ^33^ 

^2  =  '^2 

H33=^(Q,1^Q2,  ^Q33)^fP33-}(Pll^Pp2)- 


2.4  Finite  Element  Model 

This  section  describes  the  finite  element  implementation  of  the 
nonlinear  viscoelastic  constitutive  law  presented  in  Sections  2.2  and 
2.3.  Since  viscoelastic  materials  often  undergo  large  displacements 
especially  when  subjected  to  creep  type  of  loading,  the  geometrically 
nonlinear  analysis  described  in  ll,2|  has  been  incorporated  into  the 
viscoelastic  formulation. 

Invoking  the  principle  of  virtual  work  and  following  the  procedure 
outlined  in  |1,2|,  one  has 
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.  .e  sf.e.  Ady 
1  ijrs  1  rs  H  ij' 


I'ij 


)dV 


1,.  ,.(ie,j)dV  . ; 


.6u.dv  +  J  k.6u.ds, 

1  1  ^  1  1 


(2.47) 


where,  are  the  components  of  the  viscoelasticity  constitutive 

tensor.  The  rest  of  the  quantities  and  the  superscripts  in  Eq.  (2.47) 
have  the  same  definitions  as  given  in  (21.  Let  each  displacement 
increment  at  any  time  t  be  approximated  as, 


AU 


Substitution  of  Eq.  (2.48)  into  Eq.  (2.47)  gives 


((K^l  +  (K^D-au) 


■f'-\  .  'f^] 


(2.48) 


(2.49) 


where 

(K^l  =  h*  lB'-|^(M!lB^|dA,  h  =  thickness  (2.50) 

■^1 

and  (Ml  is  the  4x4  viscoelastic  constitutive  matrix  defined  in  Eq. 

(2.41).  Note  that  for  the  nonlinear  viscoelastic  case,  the  'linear' 
stiffness  matrix  1k'-|  is  not  really  linear,  but  has  nonlinearities 
imbedded  in  it  due  to  the  presence  of  the  material  kernel  functions  (g^, 
9^)  1'^  matrix  1M|.  The  nonlinear  stiffness  matrix  iK’l  is  the 
direct  result  of  the  geometrically  nonlinear  formulation  and  is  given  by 

|K^"1  =  h'  lB'’]^[T)lB'’)dA.  (2.51) 

The  definitions  of  matrices  IB*-),  (F*"},  {F“}  and  [t]  are  the  same 

as  in  (21 .  The  Cauchy  stress  components  are  computed  by  using  the 
viscoelastic  relation, 

{a|  =  (M|({e.  -  Hj)  (2.52) 

which  has  been  derived  in  Section  2.3.  For  a  geometrically  nonlinear 
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analysis,  the  vector  {e[  contains  components  of  the  Almansi  strain 
tensor 


,  3U,  3U  3U  3U 

e  =  1  f— li  +  -I _ Q  _il\  .  ,  9 

km  2  '■ax,,  ax,  ax.  3x„'  km  ‘ 

m  k  km 


(2.53) 


It  is  evident  that  Eq.  (2.49)  contains  two  possible  sources  of 
nonlinearities:  material  nonlinearity  due  to  Schapery's  law  and, 
geometric  nonlinearity  arising  from  the  large  displacement  (and  small 
strain)  formulation.  In  order  to  obtain  a  solution  to  this  nonlinear 
equation  at  any  time  step,  the  Newton-Raphson  iterative  technique  is 
used.  The  incremental  displacement  {au}^  obtained  at  the  end  of  the  r^^ 
iteration  is  used  to  update  the  total  displacement  for  the  n^*^  time 


step. 


n 

+  A  U 


(2.54) 


The  iterative  procedure  continues  until  a  convergence  criterion  is 
satisfied.  After  that,  the  solution  proceeds  to  the  next  time  step. 
Note  that  for  the  first  time-step,  the  initial  cordHion  *^0^  ^  =  0  is 
used. 


3.  MOISTURE  DIFFUSION  AND  DELAYED  FAILURE 
3.1  Governing  Equations  for  Diffusion 

The  nonlinear  Fickean  two-dimensional,  diffusion  model  presented  in 
the  present  study  is  the  one  investigated  by  Lefebvre,  et  al.  1121.  The 
diffusion  model  can  also  be  applied  for  penetrants  other  than  moisture, 
as  long  as  the  diffusion  process  is  Fickean. 

Fick's  law  for  the  two-dimensiona''  diffusion  of  a  penetrant  within 
an  isotropic  material  is  given  by 
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(3.1) 


3_  (0  1^)  +  ^  (D  1^) 
3X  '  3x'  ay  '  3y' 


where,  c  is  the  penetrant  concentration,  which  is  a  function  of  position 
and  time,  and  D  is  the  diffusion  coefficient. 

In  order  to  model  the  transport  phenomena  in  polymeric  materials, 
Lefebvre  et  al.  [121  derived  a  nonlinear  diffusion  coefficient  based  on 


the  concept  of  free  volume.  According  to  this  theory,  the  diffusion 
coefficient  for  a  polymeric  material  above  its  glass  transition 
temperature  is  given  by. 


^  -  l/f„l 


(3.2) 


where  0  is  the  diffusion  coefficient,  T  is  the  temperature,  f  is  the 
free  volume  fraction,  and  B  is  a  material  constant.  The  subscript  'o' 
denotes  values  at  the  reference  state.  It  is  then  postulated  that  the 


change  in  fractional  free  volume  is  due  to  changes  in  temperature, 
penetrant  concentration,  and  the  transient  component  of  the  mechanically 
induced  dilatational  strain.  It  is  further  assumed  that  these  changes 


are  additive,  which  is  similar  to  the  assumptions  made  by  Knauss  and 


Emri  [151.  Thus 

f  =  f  +  3c<aT  +  ae, ,  +  3-rC^.  (3.3) 

0  kk 

In  Eg.  (3.3),  a  is  the  linear  coefficient  of  thermal  expansion,  y  is  the 
linear  coefficient  of  expansion  due  to  moisture,  N  is  an  exponent  for 
the  saturated  state,  and  Aej^j^  is  the  transient  component  of  the 
mechanically  induced  dilatational  strain.  The  dilatational  strain  due 
to  the  ambient  stress  state  can  be  written  as 

®kk  "  ®kk^°^  ^  *®kk’ 

where  e|(|((0)  and  Aej^j^  are  the  instantaneous  and  transient  components  of 
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Substitution  of  Eq.  (3.7)  into  Eq.  (3.3)  and  the  result  into  Eq.  (3.2) 


yields 


D  =  T  exp  ( 


B_  I  +  yc  )  -^  (e^^^^  -  ^ 

^0  fo  ^  ^  (^kk  -  i  Vkk)  '' 


(3.8) 


From  the  viscoelastic  formulation  presented  in  Section  2,  it  is  evident 


that  the  dilatational  strain  e|^|^  is  dependent  on  the  stress  history, 
temperature  and  penetrant  concentration, 

^kk  "  ^kk^^kk’^’'^^ 

Hence,  the  two  sources  of  nonlinearity  in  Eq.  (3.1)  are  moisture 
concentration  c  and  dilatational  strain  e|^|^.  Consequently,  in  order  to 
accurately  model  the  penetrant  transport  phenomena,  the  diffusion 
boundary-value  problem  needs  to  be  solved  in  conjunction  with  the 
nonlinear  v i scoel a' c ic i ty  boundary-val ue  problem  by  using  an  iterative 
procedure.  The  same  solution  procedure  can  also  be  applied  for 
diffusion  in  polymeric  materials  where  the  plasticizing  effect  of  the 
penetrant  may  cause  the  viscoelastic  time-scale  shift  factor  to  be 


concentrat ion  dependent, 

a^^^  =  a(:,T,c)  (3.10) 
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:ne  ej-anple  of  such  a  shift  factor  definition  can  be  found  in  the  work 
of  Knauss  and  Emri  (15!,  where  the  authors  have  used  the  concept  of  free 
volume  to  define  a  shift  factor 


P  oaT  +  yC  +  xe,  , 
a(T,c,s,„)  -  exp(-  ^ 


(3.11) 


where  x  is  the  coefficient  of  the  dilatation  term.  Note  that  the 


coefficients  a  and  y  in  Eqs.  (3.8)  and  (3.11)  are,  in  general,  functions 
of  T,  c  and  ej^j^  but  have  been  assumed  to  be  constant  for  the  sake  of 
simplicity.  This  assumption  is  valid  for  temperatures  above  the  glass 
transition  temperature  and  below  the  boiling  point  of  the  penetrant. 


2.2  Fir.ite  Element  Formulation 

Pick's  law  for  two  dimensional  diffusion  in  a  homogeneous  isotropic 
material  is  given  by. 


^  (0  I£)  .  1-  (i£) 

)X  '  jx'  iy  '•iy' 


(3.12) 


subject  to  the  boundary  conditions. 


D  —  n  +  D~n  +  q  =  0onr,  ,  t:jO 
3x  X  3y  y  1 

c  =  c  on  ,  t  ?  0 


(3.13) 

(3.14) 


and  the  initial  condition, 

c=cine,t=0  (3.15) 

o 


where,  is  the  two-dimensional  region  in  which  diffusion  occurs, 
and  r  is  the  boundary  of  n. 

The  finite-element  formulation  of  Eqs.  (5. 12)-(3. 15)  can  be 
developed  following  the  variational  procedure  used  by  Reddy  (171.  The 
time  dependent  moisture  concentration  is  approximated  over  an  element 
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(3.16) 


n®  as. 


j=l 


where  N  is  the  number  of  nodes  per  element.  The  resulting  finite- 


element  equations  cast  in  a  matrix  form  are  given  by 


whene 


(3.17) 


=  J  4-^’^^dxdy 

e  e  e  e 

,  .  Sijj  3i;i .  ani.  3ii)  . 

KIV  =  J  D(-^  ^  +  —I  .:^)dxdy 

ij  ■'  „  '■ax  ax  ay  ay  ^ 


(3.18a) 


(3.18b) 


=  -  f  ^^qds  . 

.e 


(3.19) 


The  time  derivative  'Cj  is  approximated  using  the  e-family  of 


approximation. 


9{c}n+i  ^  (1  -  9)(ci 


'C  1  -  C 
'n+1  ^  'n 


n  At. 


for  0  <  9  <  1 


(3.20) 


wherefikt„  is  the  time  step  and  ci  denotes  the  value  of  Cf  at  time 
n  '  n 

t^.  Using  the  approximations  (3.22)  in  (3.18)  for  time  tj^  and 


gives 


-  ib<">i!cL  -  =  0 


(3.21) 


where 


-  (1  -  e)At^^^[K(^)l 

=  At  ,(e{F(^)}„  .  +  (1  -  9){f(")}J.  (3.22) 
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Recognizing  that  a  source  of  nonlinearity  in  the  form  of  the  diffusion 

coefficient  0  is  imbedded  in  the  matrix  the  Newton-Raphson 

technique  is  employed  to  solve  for  the  concentration  (cl  ,  at  each  time 

n+l 

step.  Note  that  for  n  =  0,  the  vector  cl^  in  Eq.  (3.21)  is  known  from 
the  initial  conditions. 

3.3  Delayed  Failure:  Uniaxial  Formulation 

When  a  viscoelastic  material  undergoes  deformation,  only  a  part  of 
the  total  deformation  energy  is  stored,  while  the  rest  of  the  energy  is 
dissipated.  This  behavior  is  unlike  elastic  material  where  all  the 
energy  of  deformation  is  stored  as  strain  energy.  Reiner  and  Weisenberg 
[51  postulated  that  it  is  this  time-dependent  energy  storage  capacity 
that  is  responsible  for  the  transition  from  viscoelastic  response  to 
yield  in  ductile  materials  or  fracture  'n  brittle  ones.  According  to 
this  theory,  failure  occurs  when  the  stored  deviatoric  strain  energy  per 
unit  volume  in  a  body  reaches  a  certain  maximum  value  called  the 
resilience,  which  is  a  material  property.  Note  that  when  there  is  no 
dissipation,  that  is,  when  the  material  is  elastic,  then  Reiner- 
Weisenberg  criterion  becomes  identical  to  the  von  Mises  criterion. 

Consider  the  single  Kelvin  element  shown  in  Fig.  1,  subject  to  the 
uniaxial  tensile  load  o(t).  The  total  strain  response  e(t)  due  to  the 
applied  stress  can  be  divided  into  two  components:  the  instantaneous 


respone  e^,  and  the 

transient  response  ej^(t), 

e(t)  =  e^  +  e^(t) 

(3.23) 

For  uniaxial  creep, 

the  applied  stress  o(t)  is  given  as 

o(t)  =  a^H(t) 

(3.24) 

where  H(t)  is  the  unit  step  function. 
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Substituting  Eq.  (3.?4)  into  Schapery's  nonlinear  uniaxial  single 


integral  law  given  by  Eq.(2.1)  and  expressing  the  transient  creep 
compliance  as 

O^(^)  =  0^  -  D^(l  -  e  (3.25) 

leads  to, 

=  %%%  ^  -  e  (3.26) 

where  li)  is  the  reduced  time  defined  in  Eq.  (2.2). 

Comparing  Eq.  (3.26)  with  (3.23),  it  follows  that 

e  =  g  D  o 
0  ^0  0  0 

-X  1 

e^(i)  =  gjD^d  -  e  )g^^^  .  (3.27) 


For  a  given  applied  stress  the  stress  developed  in  the  nonlinear 
spring  with  compliance  D^g^  is  7^  and  the  corresponding  strain 
is  e\  For  the  spring  with  the  nonlinear  compliance  D^g]^  the  stress  is 
given  by, 


s 

'l 


(3.28) 


where  the  superscript  s'  denotes  quantities  related  to  the  spring. 
From  Fig.  1  it  is  evident  that  and  e^  are  equivalent.  Hence, 


'1 


-  e 


) 


(3.29) 


The  total  energy,  W^,  stored  in  the  two  springs  over  time  t  is  (see  Hiel 
et  al .  [ 14| ) , 


^0  _t 

W  =  r  '3  de  +  a^e^dt 

0  °  'o  ^  ^ 


^i^oVo  -  e  ^  l^g^c^) 


(3.30) 
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For  a  viscoelastic  material  represented  by  multiple  Kelvin  elements 
in  series,  Eq.  (3.30)  takes  the  form. 


?  .  1  -  2 


•V\?, 


(3.31) 


According  to  the  Reiner-Weisenberg  hypothesis,  failure  occurs  when  the 
stored  energy  reaches  the  resilience  of  the  material.  Denoting  the 
resilience  as  R,  the  expression  for  the  time  dependent  failure  stress 
obtained  from  Eq.  (3.31)  for  uniaxial  stress  state  is. 


[  .  .  *Tt  ~vXpl/2 

i  %°o  "  i  ^192  I  (0,(1  -  e  )  I- 

r=l 

3.4  Delayed  Failure:  Multiaxial  Formulation 

If  Op  '2  '’3  1-1’®  principal  stresses  at  any  point  in  a 

viscoelastic  material,  then  by  definition,  the  shear  stresses  are  zero 
on  the  principal  planes.  In  order  to  simplify  the  derivation,  let  ’t  be 
assumed  that  the  viscoelastic  material  is  represented  by  means  of  a 
single  Kelvin  element  (see  Fig.  1)  in  each  principal  direction.  The 
applied  multiaxial  creep  stresses  in  the  material  principal  directions 


(3.32) 


are  given  by. 


“22  ^ 

^33  .  «3H(t). 


(3.33) 


Substitution  of  Eqs.  (3.33)  in  Eqs.  (2.31)  and  (2.32)  results  in  the 
following  expressions  for  the  corresponding  viscoelastic  strains. 


eil(t)  --  ■  0  ^‘’2  '  D  ^  '3' 

0  0 


- 

0^1(1  '  e  )g2-^ 


'  ■li  Ji  'r* 

>  1(1  -  e  M  +  fr  (1  -  e  "  )  g,c,  +  1(1  -  e 
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"oj  =  -"(^1  ^  ’2  "  ’3^  " 

- '  V 


1  U 


0  .  =  -\J 

fj 


9p(l  -  e  )('^j  ■*■  'Jp  *  o^)  +  (1  +  v)gp(l  -  e 


U 
(3.44) 


4.  NUMERICAL  RESULTS 

In  this  section  results  of  a  number  of  linear  elastic,  linear 
viscoelastic  and  nonlinear  viscoelastic  analyses  are  discussed  in  light 
of  available  experimental  or  analytical  results.  All  results  are 
obtained  using  NOVA  on  an  IBM  3090  computer  in  double  precision 
ar i thmat ic . 


4.1  Elastic  Analysis  of  a  Conposite  Single  Lap  Joint 
Renton  and  Vinson  |18|  used  a  closed  form  elastic  solution  to 
conduct  a  parametric  study  of  the  effect  of  adherend  properties  on  the 
peak  stresses  within  the  adhesive  in  a  composite  single  lap  joint.  A 
similar  pa-'ametric  study  was  carried  out  using  the  finite  element 
program  NOVA,  '"re  gecmetry,  finite  element  d  i  scret  i  zat  i  on  and  boundary 
condit'ons  for  tre  compcs'te  lap  joint  are  shown  in  Fig.  2.  For  the 
sake  of  s’^n'icity,  on'y  identical  adherends  are  considered.  Each 
adherend  is  made  up  of  seven  laminas  of  egual  thickness.  The 
orthotropic  material  properties  for  a  lamina  are  g-ven  in  Table  1.  In 
order  to  maintain  material  symmetry  about  the  laminate  mid-plane  and 
thus  eliminate  bending-stretching  coupling,  a  6' /O  /-e '/O  V-^'-VO>-/90  ply 
orientation  was  selected  for  the  analysis.  Note  that  this  type  of  ply 
orientation  places  the  e ply  immediately  adjacent  to  the  adhesive 
layer.  The  adhesive  used  is  FM-73  and  its  isotropic  linear  elastic 
properties  are  listed  in  Table  2.  The  adhesive  layer  is  modeled  using 
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Tdble  1.  Orthotropic  MdtP''idl  Froperites  for  Composite  Adherend. 

Qjj  =  46.885x10^  MPa 
OjP  -  0[3  =  4.137x10^  MPa 
^2?  '  '^33  '  14.96?xl0^  MPa 
'■21  '^3?  ■■ 

(jjj  ■  3.44/xlG^  MPa 
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Table  2. 


Isotropic  Linear  Elastic  Properties  for  FM-73. 

E  -  2.78x10^  MPa 
G  --  1.01x10^  MPa 
-  0.38 
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sixteen  eight-noded  quadr i lateral  elements  along  us  length  and  two 
elements  through  its  thickness.  A  series  of  elastic  finite  element 
analyses  is  performed  to  study  the  effect  of  ply  orientation,  lamina 
primary  modulus  dod  geometric  nonlinearity  on  the  peak  stresses 

in  the  adhesive. 

In  order  to  study  the  influence  of  ply  orientation  on  the  adhesive 
stress  distribution,  stress  analyses  were  performed  for  9  =  0°,  e  =  15’, 
e  =  45’,  and  e  =  90  respectively.  The  results  are  shown  in  Figs.  3 
and  4.  The  plots  show  the  variation  of  stresses  along  the  upper 
bondline  of  the  overlap.  The  parameter  x/c  is  the  normalized  distance 
from  the  bond  centerline  such  that  the  value  x/c  =  -1  corresponds  to  the 
left-hand  free  edge  of  the  bond  overlap.  It  is  evident  from  these 
figures  that  an  increase  in  the  ply  orientation  angle  9,  causes  the  peak 
stresses  to  increase  near  the  free  edge  cf  the  bond  overlap.  The 
adherend  with  a  0V90  '  ply  orientation  (cross-ply)  shows  a  28%  increase 
u  peel  stress  and  a  17<  inc'-ease  in  shear  stress  over  the  corresponding 
values  for  a  0  (unidirectional)  ply  orientation.  This  is  not 
surprising  since  a  cross-ply  adherend  has  a  lower  bending  stiffness 
which  results  in  a  larger  lateral  deflection  causing  higher  stress 
concentrations  at  the  overlap  ends. 

The  influence  of  the  lamina  primary  modulus  (Q^j)  on  adhesive  peel 
and  shear  stresses  can  be  seen  in  Figs.  5  and  6  respectively.  A 
unidirectional  (0-degree)  adherend  ply  orientation  is  used  for  this 
analysis.  The  two  figures  show  a  significant  increase  in  the  peak 
adhesive  stress  as  the  value  of  Qjj  decreases.  This  is  understandable 
as  a  more  flexible  adherend  would  undergo  larger  bending  and  hence 
produce  higher  stress  concentrations  at  the  overlap  ends. 
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Flfi  -  3  Innucnce  of  Ply  Orlwtatlon  on  Adheslre  Peel  Strew. 


34 


asrAN%x/c 


Figure  4  Influencx  of  Ply  Orientatioa  on  Adhesive  Shear  Stress. 
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Harris  and  Adams  [19]  conducted  large  displacement  finite  element 
analyses  on  a  single  lap  joint  with  aluminum  adherends  and  observed 
significant  reduction  in  peak  stresses  at  the  edge  of  the  adhesive  as 
compared  to  linear  results.  In  order  to  observe  the  effect  of  geometric 
nonlinearity  on  a  single  lap  joint  with  laminated  composite  adherends,  a 
large  displacement  analysis  was  performed  using  the  program  NOVA.  Due 
to  its  greater  susceptibility  to  bending,  cross-ply  laminated  adherends 
were  used  for  this  analysis.  The  results  can  be  seen  in  Figs.  7  and 
8.  The  geometrically  nonlinear  analysis  results  in  a  30%  reduction  in 
the  peak  peel  stress  and  a  15%  reduction  in  the  peak  shear  stress.  The 
horizontal  shifting  of  the  nonlinear  curves  is  due  to  the  configuration 
coordinate  update  required  by  the  large  displacement  analysis. 

4.2  Nonlinear  Viscoelastic  Analysis  of  a  Composite  Single 
Lap  Joint 

A  nonlinear  viscoelastic  analysis  of  a  lap  joint  made  of  composite 
material  was  carried  out  over  a  time  period  of  forty  hours  using  NOVA. 
The  specimen  geometry  and  the  finite  element  discrpuization  are  the  same 
as  for  the  elastic  analysis  as  shown  in  Fig.  2.  However,  instead  of  a 
uniform  end  traction,  a  uniform  end  displacement  of  0.363  mm  is  applied 
to  the  end  of  the  joint  and  is  held  constant  with  time.  The  adherends 
are  made  of  symmetric  cross-ply  laminates  whose  properties  are  given  in 
Table  1,  wnile  the  adhesive  used  is  FM-73  and  its  creep  compliance  and 
Schapery  parameters  can  be  found  in  "^able  3. 

Figures  9  and  10  show  the  variation  of  shear  stress  and  shear 
strain  respectively  across  the  entire  bond  length  over  a  period  of 
hours.  The  sharp  peak  on  the  left-hand  edge  is  due  to  the  presence  of  a 
re-entrant  corner  and  also  due  to  the  difference  in  material  properties. 
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Figure  7 .  Influence  of  Geometric  Nonlinearity  on  Adhesive  Ped  Stress. 
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Figure  8.  ’nduence  of  Geometric  Nonlinearity  on  Atlliesive  Shear  Stress. 


Table  3.  Material  Data  for  FM-73  Unscrimmed  at  30"C. 


Elastic  Compliance,  D^;  360  x  lO'^/MPa 

Poisson's  Ratio,  v:  0.38 

Coefficient  of  Thermal  Expansion,  a:  6.6  x  iO'^  m/m/“K 
Prony  Series  Coefficients: 

0^ 

O3 

D4 
O5 


=  11.05xl0'®/MPa 

=  10  secs. 

=  12.27xl0-^/MPa 

=  10^  secs. 

=  17.35xl0'^/MPa 

=  10^  secs. 

=  21.63xl0'^/MPa 

=  10^  secs. 

=  13. 13xI0-^/MPa 

=  10^  secs. 

=  41.78xl0'^/MPa 

■g  =  10^  secs. 
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Figures  11  and  12  provide  a  ciose-up  view  of  the  sheir  stress  and  st'-im 
gradients  at  the  free  edge.  As  might  be  expected,  the  shear  stress 
undergoes  relaxation  v<hich  results  in  a  26%  decrease  in  tne  peak  value 
at  the  left  hand  edge.  The  stresses  have  been  normalized  with  respect 
to  an  average  shear  stress  value  of  4.6  N/mm^.  The  peak  shear  strain, 
however,  shows  an  increase  of  35%  over  the  same  period  of  time. 

Similarly,  Figs.  13  to  16  revea’  that  while  the  peak  values  of  the  pe? 
and  axial  stresses  decrease  by  26%  and  32%  respectively,  the 
corresponding  strains  show  a  respective  increase  of  63%  and  6%.  The 
reason  that  the  strains  increase  with  time  even  though  the  joint  end 
deflection  remains  fixed,  is  because  the  adherends  are  modeled  as 
elastic  continuums.  As  the  stresses  in  the  adhesive  relax  with  time, 
the  elastic  adhererds  deform  to  attain  a  new  equilibrium  configuration 
and  this  leads  to  an  altered  state  of  strain  wUhin  the  adhesive. 

Hence,  it  is  very  important  that  the  elastic  nature  of  the  adhererds  ne 
taken  into  account  in  an  analysis.  Also,  the  significant  increase  ’n 
adhesive  strains  with  time  is  a  viscoelastic  pheno'.Lnrn  and  therefcre  t 
cannot  be  precicteo  by  means  of  a  purely  elastic  analysis,  "ris  fact 
emphasizes  the  need  to  model  the  adhesive  ayer  as  a  viscoelastic  medium 
in  order  to  be  able  to  predict  the  long-term  durability  of  a  bended 
joint. 

4.3  Nonlinear  F ickean  Diffusion  in  Polystyrene 

In  order  to  validate  the  diffusion  model  implemented  in  NOVA  and 
discussed  in  Section  3,  results  from  a  nonlinear  diffusion  analysis 
presented  in  (201  are  used.  The  test  problem  involves  unsteady  sorption 
of  a  penetrant  in  a  semi- inf inite  medium  for  a  diffusion  coefficient 
that  is  an  exponential  function  of  penetrant  concentration,  that  is. 
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Figure  i  3  Variation  of  Peel  Stress  with  Time  Near  the  Free  E<ige. 
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J  ^  e^p  (kC/CQ).  Finite  element  predictions  -vere  obtained  for  k  = 
F.olA  and  k  =  3.912  and  the  results  were  compared  with  the  pubMshed 
results  represented  by  the  solid  lines  in  Fig.  17.  Excellent  agreement 
is  observed  for  the  two  values  of  the  coefficient  k. 

Levita  and  Smith  |21]  conducted  experiments  to  study  gas  transport 
in  polystyrene  S' d  round  that  the  diffusion  coefficients  for  gases 
decreased  with  time  when  the  polystyrene  film  was  subject  to  a  constant 
uniaxial  strain.  This  effect  was  attributed  to  the  continuous  free 
volume  recovery  (densi f ication)  in  the  polystyrene  specimen  at  constant 
strain.  The  study  also  indicated  that  larger  free  volume  elements 
Increase  in  si^e  faster  than  the  smaller  ones  as  volume  recovery 
0'-ogr.:>sses .  Using  the  results  published  in  [211  as  a  guideline,  NOVA 
•■.as  uivo  to  stud;,  the  time  dependence  of  the  diffusion  coefficient  for 
;  I'to'-- 1 -C' ice  gas  -n  a  polystyrene  film  at  constant  strain.  For  this 
tn?  '--mc^dture  and  moisture  concentration  effects  presented  in  Eq. 
].“)  Aprc-  '■pc’ected,  resulting  in  a  dit^fusion  coefficient  that  is 

a  fu'-'  tier  of  the  transient  component  of  the  dilatational  strain 
t,.rr,  -s  a  measure  of  the  change  in  the  free  volume.  Figure 
S'":.'.s  the  var-:dt'on  of  the  diffusion  coefficient  with  time  for  three 
;  strain  levels.  The  material  properties  for  polystyrene  which 

<vere  uttained  from  |27l  are  given  in  Table  4.  From  Fig. 18  it  is 
evident  that,  independent  of  the  strain  level,  the  diffusion  coefficient 
^eauh^-'S  a  peak  value  at  around  t  =  1  hour  and  then  slowly  decays  to  the 
rpfor-Tree  value,  This  behavior  can  be  attributed  to  an  initial 

’nrrease  in  free  volume  due  to  the  application  of  the  uniaxial  strain, 
^}llo.ved  hy  a  continuous  recovery  in  free  volume  (dens i f i cat ^ on)  at  a 
o^std'-t  st^'ain  as  the  polystyrene  film  undergoes  relaxation.  A  larger 
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Table  4.  Material  Properties  for  Polystyrene  at  50  C. 


Bulk  Compliance: 

Mq  =  l.PxlO-'^/MPa 
M^  =  0.2896xlO-^/MPa 
M2  =  0.2246xl0-^/MPa 
M3  =  0.3721xl0-^/MPa 
M^  =  0. 1354xlO"^/MPa 
Shear  Compliance: 

Jq  =  l.OxlO'VwPa 
=  2.16/MPa 
J2  =  2.92/MPd 
J3  =  1.38/MPa 
J4  =  2.88/MPa 
J5  =  2.31/MPa 
Jg  =  3. 59 /MPa 
J7  =  0.648/MPa 


=  1.515x10^  sec. 
=  1.515x10^  sec. 
T3  =  1,515x10'^  sec. 
=  1.515x10^  sec. 


=  1 

.515x109 

sec. 

=  1, 

.515xl0l° 

sec. 

=  1, 

.515x10^^ 

sec. 

=  1, 

.515x10^^ 

sec. 

=  1. 

,515x10^'^ 

sec . 

=  1. 

.5i:xl0^'’ 

sec. 

=  1. 

515xlo'-^ 

sec . 

Reference  free  volume  f^  =  0.033 
Oiffusion  coefficient  ^  9x10'^  mm^/sec 
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applied  strain  produces  larger  initial  dilatation,  and  this  results  in  a 
higher  peak  in  the  diffusion  coefficient.  Figure  18  also  reveals  that 
the  time  rate  of  free  volume  recovery,  and  hence  the  time  rate  of 
decrease  in  the  diffusion  coefficient,  is  proportional  to  the  applied 
strain  level. 

The  influence  of  penetrant  molecule  size  on  the  diffusion 
coefficient  for  gases  in  polystyrene  was  studied  by  varying  the 
magnitude  of  the  material  parameter  B  in  Eq.  (3.8).  The  temperature  and 
strain  were  held  constant  at  50' C  and  1.8%  respectively.  The  prediction 
obtained  from  NOVA  are  shown  in  Fig.  19  for  two  values  of  B.  The  faster 
rate  of  decrease  in  the  diffusion  coefficient  for  a  higher  value  of  B 
implies  that  the  larger  free-volume  elements  decrease  in  size  faster 
than  the  smaller  ones  as  volume  recovery  progresses.  The  NOVA 
predictions  are  qualitatively  in  good  agreement  with  the  results 
presented  in  (2I|. 

When  a  polymeric  material  is  in  the  rubbery  state,  equilibrium  is 
reached  very  rapidly  in  response  to  variations  in  temperature,  stress 
and  penetrant  conentration.  By  contrast,  a  material  in  the  glassy  state 
is  not  in  thermodynamic  equilibrium  and  tne  response  of  the  free  volume 
to  changes  in  external  conditions  is  delayed.  This  metastable  state 
causes  the  free  volume  to  slowly  collapse  with  time  until  equilibrium  is 
reached.  This  phenomenon  is  known  as  physical  aging  and  causes 
relaxation  processes  to  take  place  over  a  longer  time.  Struik  [231 
proposed  that  for  a  material  in  the  glassy  state,  effective  time  x.  is 
related  to  actual  time  t  by, 

t  t  t. 

0  ^^e  ^ 
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where  is  the  aging  time  at  the  start  of  service  life  or  testing 
and  u  is  a  constant  such  that  0  <  u  <  1.  For  such  a  material,  the 
definition  of  reduced  time  given  by  Eg.  (2.2)  is  no  longer  valid  and 
should  be  modified  to. 


where  a^,  is  the  shift  factor, 
a  I 

The  effect  of  physical  aging  on  the  diffusion  coefficient  for 
carbon-dioxide  gas  in  polystyrene  was  studied  by  implementing  Eqs.  (4.1) 
and  (4.2)  in  NOVA.  The  values  of  temperature,  strain  and  t^  were  set  at 
50‘'C,  l.Q%  and  24  hours  respectively.  Figure  20  shows  that  an  increased 
physical  aging  denoted  by  a  higher  value  of  the  parameter  u,  causes  the 
diffusion  coefficient  to  decay  slower  than  the  one  for  which  •„  is 
lower.  This  behavior  is  expected  since  increased  physical  aging  causes 
the  free  volume  recovery  to  take  place  over  a  longer  period  of  time. 

Note  that  when  there  is  no  physical  aging,  .  and  t^  are  equal  to  zero 
and  \  is  identically  equal  to  t. 

4.4  Sonl  iP.ear  Viscoelaszic  .Analysis  of  a  Butt  Joint 
Including  Moisture  Diffusion 

The  effect  of  a  change  in  the  free  volume  of  a  polymer  on  its 
viscoelastic  response  was  discussed  by  Knauss  and  Emri  [151.  They  used 
the  unifying  concept  of  the  free  volume  by  considering  that  fractional 
free  volume  depends  on  three  variables;  temperature  T,  moisture 
concentration  c,  and  mechanically  induced  dilatation  9.  Lefevbre  et  al. 
[121  extended  the  free  volume  concept  to  define  a  nonlinear  diffusion 
coefficient,  which  results  in  a  coupling  between  the  viscoelasticity  and 
the  diffusion  boundary  value  problems  (see  Section  3.1).  The  influence 
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Figure  2'  EITect  of  Physical  Aging  on  the  Diffusion  CocfTicient  for  Polystyrene. 


of  th's  coupling  on  the  viscoelastic  response  anc  no^sture  ii^'fuSion 
within  the  adhesive  layer  of  a  butt  joint  was  investigated  by  using  the 
program  NOVA.  The  specimen  geometry  and  finite  element  discretization 
are  shown  in  Figs.  21  and  22,  respectively.  A  uniform  axial 
displacement  of  0.1  mm  is  applied  at  the  end  of  the  joint  and  is  held 
constant  with  time.  The  adherends  are  made  of  aluminum  and  the  adhesive 
used  is  polystyrene.  The  various  material  properties  are  listed  in 
Tables  4  and  5.  The  selection  of  polystyrene  as  an  adhesive  was 
prompted  by  the  fact  that  it  is  one  of  the  few  polymeric  materials  that 
have  their  viscoelastic  properties  and  diffusion  parameters  adequately 
documented.  The  normalized  moisture  concentration  at  the  free  edge  of 
the  aohesive  layer  is  unity,  and  the  initial  concentration  th'^oughout 
the  adhesive  layer  is  zero.  The  tests  are  conducted  at  the  reference 
temperature  of  50“C. 

Figure  23  shows  the  moisture  concentration  profi’es  within  the 
adhesive  layer  at  three  different  times  when  there  is  no  coupling.  In 
this  case  the  diffusion  coeffic'ent  remains  constant  with  time,  that  is, 
0  =  Figui'r)  24  sho.-.s  the  mo  store  concentration  profiles  for  the 

case  where  there  is  viscoelastic  ccuoling  only,  that  's,  when  the 
diffusion  coefficient  depends  on  tne  transient  component  of  tne 
dilatational  strain.  Figure  25  depicts  the  case  where  there  is  full 
coupling,  that  is,  the  aiffusion  coefficient  is  a  function  of  the 
dilatational  strain  and  the  moisture  concentration  at  any  given  point  in 
the  adhesive.  Conversely,  the  viscoelastic  shift  factor  is  now  a 
function  of  the  dilatational  strain  and  the  moisture  concentration. 
Figure  26  presents  the  results  for  each  of  these  three  cases  for 
comparison  at  time  t  =  8  hours.  From  these  figures  it  is  evident  that 
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Figure  24  Moisture  Pronics  Within  the  Adhesive  For  Only  Viscoelastic  Coupling. 
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Figure  26  Influence  of  Coupling  on  Moisture  Profilea  at  Time  "  8  Hours. 
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Table  5.  Properties  for  Elastic  Analysis  of  a  Butt  Joint. 


Materials 

E(MPa) 

Steel 

2.07x10^ 

0.29 

Aluminum 

0.7x10^ 

0.33 

Eponal 

5.8x10^ 

0.33 

Rigid  Epoxy 

2.2x10^ 

0.33 
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the  effect  of  coupling  is  to  accelerate  moisture  diffusion  in  the 
adhesive.  The  mechanically  induced  dilatation  together  v*itn  the 
swelling  due  to  moisture  sorption  results  in  a  higher  free  volume 
fraction  within  the  adhesive  which,  according  to  Eg.  (3.9),  causes 
diffusion  to  proceed  faster  over  the  same  period  of  time.  It  is  to  be 
noted  that  in  Fig.  26  the  curves  become  less  concave  as  the  coupling 
increases,  which  is  in  good  agreement  with  the  results  published  in 
[20). 

Figures  27  to  30  show  the  variation  of  the  stresses  and  strains 
with  time  within  the  adhesive  layer  in  the  butt  joint  when  there  is  no 
coupling  due  to  moisture  induced  swelling.  Mathematically,  this  implies 
that  Y=  0  in  Eg.  (3.3)  and  (3.11).  From  Figs.  27  and  28  it  is  evident 
that  the  stresses  do  not  relax  significantly  over  the  time  period  of  the 
analysis.  This  is  because  the  elastic  adherend  acts  as  a  spring  causing 
the  adhesive  to  creep  even  though  the  joint  end  displacement  remains 
fixed.  However,  there  is  a  slight  relaxation  in  the  normal  stress  as 
one  moves  towards  the  center  of  the  bond.  The  large  increase  in  the 
strains,  as  sho^n  in  Figs.  29  and  30,  is  due  to  the  creep  caused  by  the 
strain  recovery  in  the  elastic  adherend.  This  observation  is  supported 
by  Fig.  31  which  shows  that  the  normal  strain  in  the  adherend 
immediately  adjacent  to  the  interface  undergoes  significant  reduction 
with  time.  The  decrease  in  the  adherend  normal  stress,  as  shown  in  Fig. 
32,  reflects  the  concurrent  stress  relaxation  that  occurs  in  the 
adhesive  and  triggers  the  strain  recovery  in  the  adherend. 

Figures  33  to  36  show  the  effect  of  moisture  induced  swelling  on 
the  viscoelastic  stresses  and  strains  in  the  adhesive  layer. 
Mathematically,  this  means  y  has  a  nonzero  value  in  Egs.  (3.8)  and 
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Figure  d  Variation  of  Normal  Stress  in  the  Adhesive  With  Time  For  Viscoelastic 
Coupling. 
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(lie  Adhesive  With  Time  For  Viscociastie 


Figure  32  Variation  of  Nonnai  Strain  in  the  Adhcrcnd  With  Time  For  Viscodaslic 
Coupling. 
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Figure  34  Variation  of  Shear  Stress  in  the  AUhcsive  With  Time  For  Vlscodastic  and 
Moisture  Coupling. 
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V ariation  of  NormaJ  Strain  in  the  Adhesive  With  Time  I'or  ^'iscoelastic  and 
Moisture  Coupling. 
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Figure  36  Vaiialton  of  Shear  Strain  In  the  Adhesive  With  Time  For  VlscoeiasUe 
Mobture  Coupiing. 


(3.11).  The  dctudl  vdlue  of  >  selected  for  this  stuc>  is  0.001.  For 
this  value  of  y,  the  moisture  absorbed  causes  large  sv^ell'nq  strains 
within  the  adhesive,  which  increase  in  magnitude  as  the  diffusion 
progresses.  This  moisture  induced  swelling  strain  causes  a  -eduction  in 
the  mechanically  induced  normal  strain  and  hence  a  icwer  value  fo'-  the 
normal  stress  in  the  adhesive.  This  effect  can  be  observed  in  Fig.  33 
where  progressive  swelling  has  caused  a  25%  reduction  in  the  peak  normal 
stress  over  a  period  of  8  hours.  It  is  interesting  to  note  that  the 
difference  between  the  two  stress  curves  diminishes  as  one  moves  towards 
the  center  of  the  bond,  ^his  behavior  is  expected  since  there  is  very 
little  moisture  near  the  center  of  the  bond  and  sc  the  stress  reduction 
is  primarily  due  to  viscoelastic  relaxation.  Tne  large  increase  in  the 
adhesive  strains,  as  seen  in  Figs.  3S  an: 36  ,  is  cue  to  the  adhererd 
acting  as  a  elastic  sor'ng. 

Figure  37  snows  thn  ''■•‘'uence  of  the  moisture  coefficient  ,,  on  the 
normal  stress  in  the  adhesive  ia>e-  after  eiqnr.  hours  of  sorotion. 
can  be  seen,  the  swelling  inpucea  for  ,  =  O.GQl  reuu'ts  in  a 
s i gn i f i cant iy  lower  normal  stress  near  tne  fray  edge  as  compared  to  tne 
case  where  .  ^  0.  Away  from  the  free  edge,  th^i'  two  stress  curves  aorear 
to  merge  as  one  moves  towards  t'-p  inte-’'-  of  tre  iiord.  This  is  oecausc 
the  low  moist.jre  concentrat  i cns  oresent  in  the  bord  interior  'S 
i  nsuf  f  i  c  ’ -..Mt  to  cause  any  significant  -eduction  in  the  normal  stress  due 
to  Swe  1  !  ' nq . 

4.5  .Of' ]  ri  q.’.i  (jf  -i  3':t  :  Joint 

The  theory  orespnted  in  Seis.  3.3  and  3.4  was  applied  to  predict 
viscoelastic  creep  failure  w'thin  the  adnesive  in  a  butt  joint.  The 
specimen  geometry  and  the  finite  element  d i scret i cat  ion  are  the  same  as 
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Adhesive  NormaJ  Stress  at  Time 


shown  in  Figs,  i-'l  and  ?2,  respectively.  ihe  adherenJ  is  made  cf 
aluminum  and  its  material  properties  are  given  in  Table  5.  ihe  adhesive 
used  is  FM-73  and  its  tensile  creep  compliance  is  listed  in  Table  3. 

The  failure  parameter  (R)  for  FM-73,  also  known  as  1*^0  rcsii'cnce,  was 
obtained  by  computing  the  area  under  the  stross-stna i n  curve  prcsc'sed 
in  [24].  This  procedure  yielded  a  value  of  the  resilience  as  1.3 
N.mm/mm^.  Note  that  the  area  under  the  visco-plastic  yield  plateau  was 
not  included  in  computing  the  value  of  R.  According  to  the  Reiner- 
Weisenberg  theory,  failure  occurs  when  the  stored  energy  per  unit  volume 
in  the  body  reaches  the  resilience  R,  for  the  material.  Using  this 
postulate  as  a  failure  criterion,  NOVA  was  utilized  to  predict  failure 
in  the  adhesive  layer  of  tne  bjtt  joint  suoject  to  a  constant  uniaxial 
tension.  The  influence  o'  applied  stress  level  on  decayed  failure  .vus 
studied  by  using  a  stress  level  of  69,  6C.  and  54  MPa  respectively.  In 
all  three  cases,  fa'ld^e  •■.as  initiated  in  the  adhesive  element  ccated 
right  at  the  free-edge  and  immediately  adjacent  ■;  tne  interface.  It 
was  also  observed  that  the  direction  of  the  plane  j-'  fr  lure  was  always 
inclined  at  an  angle  of  13  ,  counter-c lock ’ se  to  the  x-axis.  Since  the 
direction  of  "ailure  coincides  with  the  oi  ection  of  p''incipd1  stress, 
it  is  evident  that  a  multiaxiai  state  o^'  stress  exists  rear  the  free 
edge,  even  though  the  applied  stress  is  uniaxial.  Figure  3S  shows  the 
variation  of  normal  (o-  creep)  strain  with  time  at  30  C  for  the  e'enent 
in  which  failure  is  first  initiated.  The  ^ight  hand  terminat icn  point 
on  the  curves  indicate  the  noint  at  which  fa’ lure  occured.  it  is 
observed  that  for  an  applied  stress  level  of  69  MPa,  the  time  to  fa  lure 
(tp)  is  1.5  secs.  In  other  words,  for  this  stress  level,  failure  occurs 
almost  i nstanlaneous 1 V .  tor  an  applied  stress  of  60  MPa,  t-  inc'^eases 
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Figure  38  Creep  Strain  in  FM-73  for  DifTerent  Applied  Stress  Levels. 
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to  400  secs.  Reducing  the  applied  stress  to  54  MPa  results  in  a  tiae  to 
failure  of  approximately  10  hours.  These  results  are  qual itativel>  in 
good  agreement  with  the  results  presented  by  Bruller  (6!  fo--  PMMA. 

From  the  above  observations  it  is  clearly  evident  that  for 
.'iscoelastic  polymers  like  FM-73,  the  time  to  failure  depends  strongly 
on  the  applied  stress  level.  Figure  o9  shows  the  evolution  of  stored 
energy  with  time  for  different  stress  levels.  For  very  high  applied 
stress  levels,  almost  all  the  strain  energy  is  conserved  as  stored 
energy  and  failure  occurs  almost  immediately.  For  intermediate  levels 
of  applied  stress,  viscoelastic  creep  causes  a  part  of  the  strain  energy 
to  be  dissipated.  As  a  resu't,  only  a  fraction  of  the  total  strain 
energy  is  conse’^ved  as  stored  energy.  Consequently,  the  stored  energy 
builds  up  slowly,  ana'tfjCi^s  t .)  a  “leaking  vessel",  nesulf’ng  in  de'ayed 
failure.  For  st'-ess  ’ev>^'  f^at  is  O'e’ow  a  certain  threshold 

value  for  a  given  mat-'-'i',  f'o  -Csslpiteo  eneigy  nay  constitute  a  lange 
fraction  of  the  t''ta'  -'n,-'-jy.  In  that  case,  tne  sto^'et  enercv 

would  increase  toe  s'e-v. '>  to  exceed  the  res  i  1  ierc'.'  .i  *  le  material  cve'^ 
any  realistic  lengt^^  ct  t'^ne,  and  nence  tnere  ..oulo  be  no  failure  even 
if  the  applied  stress  acts  indefinitely. 

h.  summary  and  conclusions 

5.1  Ger.t^ral  Su.r.niary 

The  adhesive  layer  is  modeled  using  Schapery's  nonlinear  single 
integral  constitutive  law  for  uniaxial  and  multiaxial  stress  states. 

The  effect  of  temperature  and  stress  level  on  the  viscoelastic  response 
is  taken  into  account  by  a  nonlinear  shift  factor  definition.  Penetrant 
sorption  is  accounted  for  by  a  nonlinear  Fickean  diffusion  model  in 
which  the  diffusion  coefficient  is  dependent  on  the  temperature, 
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penetrant  concentration,  and  the  dilatacional  stram.  A  ce'ayed  failure 
criterion  based  on  the  Reiner-Wei senberg  failure  theory  oas  also  been 
implemented  in  the  finite  element  code.  The  program  is  validated  oy 
comparing  the  present  results  with  analytical  and  experimental  resuUs 
available  in  the  literature.  The  program  capaoility  has  been  extended 
to  account  for  laminated  composite  adherends  and  adnesives  with  a  time 
dependent  Poisson's  ratio.  In  general,  the  computer  program  NOVA,  is 
believed  to  provide  accurate  predicticns  over  a  wide  range  of  specimrp 
geometries,  external  loads,  and  environmental  conditions. 

5.2  Conclusions 

The  results  presented  in  Section  4  underscore  the  importance  of 
modeling  the  adhesive  in  a  bended  joint  as  a  viscoelastic  material. 

This  allows  the  analyst  to  predict  the  large  increments  in  adnesive 
strains  that  occur  with  time  and  cannot  be  predicted  by  a  pure’y  elastic 
analysis.  Furthermore,  other  events  (such  as  moisture  diffusion  an;] 
delayed  failure),  that  are  highly  relevant  for  ponded  joint  analysis, 
cannot  be  accurately  predicted  unless  viscoelasticity  is  taken  into 
account.  At  high  stress  levels,  nonlinear  viscoelast'c  effects  can 
produce  creep  strains  that  are  significantly  larger  than  the  linear 
viscoelast’c  predictions  and  such  effects,  therefore,  should  be 
accounted  for.  The  effect  of  change  in  Poisson's  ratio  with  time  in 
some  polymers  have  a  significant  bearing  on  tne  final  response  and  must 
be  taken  into  account  in  erde'"  to  obtain  accurate  results. 
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